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We study the quantum dynamics of an ultracold atomic gas in a deep optical lattice within an 
optical high-Q resonator. The atoms are coherently illuminated with the cavity resonance tuned 
to a blue vibrational sideband, so that photon scattering to the resonator mode is accompanied by 
vibrational cooling of the atoms. This system exhibits a threshold above which pairwise stimulated 
generation of a cavity photon and an atom in the lowest vibrational band dominates spontaneous 
scattering and we find a combination of optical lasing with a buildup of a macroscopic population 
in the lowest lattice band. Including output coupling of ground-state atoms and replenishing of hot 
atoms into the cavity volume leads to a coherent, quantum correlated atom-photon pair source very 
analogous to twin light beam generation in a nondegenerate optical parametric oscillator. 



Soon after the experimental realization of Bose- 
Einstein condensation of ultracold atoms, a coherent 
atomic beam was produced by continuous RF out- 
coupling of atoms from a BEC in a magnetic trap 
Although duration of such quasi continuous atom las- 
ing is limited by the total number of condensed atoms 
it constituted the first realization of an atom laser, soon 
followed by closely related implementations [2|, |3| • Even 
earlier in theory atom lasers attracted a great deal of 
interest and a variety of proposals to implement contin- 
uous atom lasers were made In practice a central 
hurdle towards a cw implementation is the presence of 
near resonant light in the cooling process, which hampers 
the buildup of a degenerate matter wave state by reab- 
sorption heating if densities sufficient for evaporation are 
involved. One way to avoid this is spatial separation of 
the laser and evaporative cooling stages in a long beam 
geometry. This should allow continuous operation [5| but 
the requirements to reach sufficient atomic flux for fast 
evaporation in the final stage are challenging. 

In this Letter we study a scheme involving concurrent 
generation of a coherent atomic beam and a laser beam 
in a single device. It consists of laser cooled atoms in a 
far off-resonant optical lattice with some spatial region 
of the lattice confined within an optical cavity. Lasing 
on an anti-Stokes vibrational Raman transition should 
dominate spontaneous emission and replace evaporation 
in the last cooling stage to reach atomic quantum de- 
generacy 0, [IH . Thus the cavity solves the reabsorption 
problem and lowers density requirements as no collisional 
thermalization is needed. 

In principle it has been noted already some time ago 
that ground state cooling of a single trapped particle can 
be enhanced by employing a cavity Q. This scheme, 
proposed to operate in the bad cavity regime, has been 
recently generalized and analyzed in more detail also for 
narrow bandwidth cavities . In both cases the cavity 
acts as a passive element. As low cavity photon numbers 
are required this cannot be simply scaled to larger vol- 
umes with many particles. Nevertheless in a transient 
regime collective enhancement of cavity side band cool- 
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FIG. 1: Schematic of the atom-photon pair laser. Hot atoms 
tunneling into the resonator accumulate in the lowest (laser) 
band while a coherent cavity field builds up. 



ing was predicted in a low excitation regime [loj]. In a 
parallel development we studied light forces in a single- 
atom laser involving an incoherently pumped pointlike 
atom in a cavity. When optical gain is tied to motional 
cooling, lasing, cooling and trapping occur together 11 1 
with final temperatures below the standard cavity cool- 
ing limit [13] ■ 

In this work we combine these ideas and study a Ra- 
man laser operated on the vibrational anti-Stokes tran- 
sition of the quantized trapped atom states of an opti- 
cal lattice. This combines gain enhanced cavity cooling 
with collective amplification of ground-state occupation 
and suppressed spontaneous emission. Above thresh- 
old we expect buildup of a coherent cavity field and a 
macroscopic atomic population utilizing Bosc enhance- 
ment from both atoms and photons. 

Note that for large detuning of the Raman transition 
from the upper atomic state, cavity loss as the main 
dissipation channel acts like a diode on the vibrational 
transition. Cavity photons and ground-state atoms are 
produced in pairs but photons have a smaller lifetime 
than atoms in the lowest band, so that a much larger 
atom population will build up. Atomic out-coupling is 
required to prevent self termination of the laser. In ad- 
dition atomic heating will enhance the photon gain and 
accelerate cooling. 
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FIG. 2: (a) Transition scheme for the Raman gain. When 
the cavity-laser frequency difference matches the energy gap, 
heating (gray path) is suppressed, (b) Illustration of the 
atomic dynamics. 



Let us now construct a simple model Hamiltonian. 
First we have the two lasing modes for the photons and 
atoms, i.e. a single resonator mode for the light and the 
lowest lattice band for the atom laser mode. Some pre- 
cooling assures sufficiently many atoms in the lowest two 
vibrational states of the lattice, which at this point is a 
ID lattice of atoms transverse to the cavity mode (see 
Fig- H]) a s discussed as well in Refs. 13|, LLJ]. In the sim- 
plest case the lattice laser also takes the role of the Ra- 
man pump and thus is assumed to be red detuned from 
the cavity mode by the energy difference between the two 
atomic bands. Raman gain induces motion-photon cou- 
pling so that in the Lamb-Dicke limit [l5| each scattering 
event induces a concurrent change of the atomic motional 
state and photon number by one. An illustration of the 
Raman gain transition can be found in Fig. [5] (a) where 
m indicates the motional quantum number of the atom. 

A corresponding effective single-particle Hamiltonian 
for a far detuned Raman transition coupling the cavity 
mode and the motional atomic dynamics reads [l6| 



H = ^- + V sin 2 (kx) 
Zm 



h(A c + U )a 1 a 



-i\/hVoUo sin(fca;) (a^ — a) 



(1) 



Here a is the operator annihilating a cavity photon while 
Uq describes the light shift introduced by the atom-field 
coupling. The depth of the lattice potential is Vq and A c 
denotes the detuning of the cavity resonance from the 
lattice field. We can now rewrite Hq using the ladder- 
operators B and to get 

H = ilB^B + (A c + U )a)a + iri{<J - a) (B + B+) , (2) 



where H = 1, fi = ky/2V /m, r\ = I^/WqUo, and I -C 1 
is the Lamb-Dicke parameter. 

When the cavity resonance is tuned to the blue side of 
the laser according to A c +C/ = ^, each scattered photon 
picks up an extra energy Ml by which the atomic motion 
is cooled. For sufficiently large f2 (deep potential) the 
offresonant heating process indicated by the gray path in 
Fig. [2] (a) can be neglected. Omitting the corresponding 
terms we can rewrite the Hamilton operator @ in the 



(3) 



Notice that H conserves the sum of photons and vibra- 
tional quanta and besides cooling also the reverse process 
where a photon is absorbed from the cavity is possible. 
However, when a photon leaves the cavity, it takes away 
energy and entropy from the system which cannot be re- 
absorbed. Hence the atom will eventually end up in the 
ground state with the mode in the vacuum state (loj . 

Let us now generalize the single-particle Hamiltonian 
((T|) into a second quantized form for many atoms follow- 
ing Ref. 17] , where we will consider the lowest two levels 



of the lattice wells corresponding to the first two levels 
(called b and c) in the Wannier expansion of the atomic 
field operator [see Fig. [5] (b)]. We assume here that the 
light shift Uo introduced by the atoms is small so that we 
can drop terms proportional to Uq in the second term of 
Eq. Thus the interaction Hamiltonian reads 
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ir\ (a)b^c — a6c^) 



(4) 



Here, b and c are the bosonic field operators for the 
atomic laser and the source mode, respectively. The term 
a^b^c describes a transition of an atom from the source 
mode to the laser mode while a photon is created in the 
cavity. Note that Eq. (0| is equivalent to nondegener- 
ate parametric amplification, which has been extensively 
studied for squeezed light generation (l8j |. For parame- 
ters yielding squeezing we thus expect nonclassical atom- 
photon correlations. Upconversion (atomic vibrational 
heating) only occurs when both modes are excited and 
light loss causes a net downward transfer of atoms. Here 
we omitted direct atom-atom interaction terms, as we 
are dealing here with an all-optical setup, where we can 
expect only small occupation numbers per site. 

An important ingredient in the model is a proper and 
consistent description of pumping (atom injection) into 
the source mode c. The simplest mechanism is atomic 
tunneling between lattice sites outside and within the 
cavity mode. As in general the tunneling in the lower 
atom-laser band b is much slower than in the source 
mode, it can be neglected so that only excited band atoms 
move. This is a coherent process and a possible Hamil- 
ton operator has the form ir^c — iTc' where the operator 
T models all coupled sites of the lattice potential out- 
side the cavity. When their number is large enough and 
they are sufficiently occupied, their quantum state can 
be approximated by a coherent state. Tracing out the 
corresponding levels from the density matrix, we get an 
effective pump Hamiltonian of the form 



i/j. (c^ — c) . 



(5) 



Within this approximation the number of atoms in the 
system (modes b and c) is not conserved. Nevertheless 
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the total atom number is conserved so that the atomic 
pump phase enters as an artifact with no observable con- 
sequences. Equation ([5]) gives thus a good approximation 
to injection by tunneling from a coherent reservoir cover- 
ing the essential physics in a self-consistent way. Similar 
to injection atoms are coupled out from the lower band. 
We will model this simply by a decay rate (analogous 
to cavity loss) which describes the coupling to an empty 
mode continuum as it would be realized by selective RF 
out-coupling in an extra magnetic field gradient. 

Figure [3] (b) sketches the atomic dynamics includ- 
ing pump, excitation transfer, and damping of the laser 
mode (dark arrows) . To be more realistic we enlarge this 
idealized dynamics by upper-state losses and incoherent 
atomic transitions where no photon is emitted into the 
resonator mode. This models residual spontaneous emis- 
sion and atomic background collisions heating the atomic 
motion. Here, affected ground-state atoms are incoher- 
ently transferred to the source mode while upper-state 
atoms are just removed from the system as indicated by 
the gray arrows in Fig. [2] (b). 

In total the master equation reads 

P = -i [-ffgain + H pump , p] + Cp, (6) 

where the Liouville operator Cp consists of the loss terms 

C k p = Kk ( [k, p + [k p, A;*] ) , k = a, 6, c (7) 

and the heating term of the atom laser mode 

C bcP = K hc ([b J,ptfc] + [bJp,tfc]). (8) 

Following standard procedures this operator equation 
can be converted into a c-number Fokker-Planck equa- 
tion for the generalized P representation, which yields 
the following set of Langevin equations, 

a = — /s a + J7/3'7 + £ ) (9a) 
P = - [Kb + n bc (7 f 7 + 1)] $ + T) a f 7 + £b , (9b) 
7 = - («c - Kbc/3 t /3) 7 - rj a/3 + p + £ c , (9c) 

where £ a , £b, and £ c are stochastic forces with zero mean. 

The semiclassical approximation to the dynamics 
amounts to ignoring the stochastic nature of Eqs. 
so that for ij 2 > K a Kb c an oscillation threshold emerges 
at 



K a (Kb + n bc ) 

Pth = Kc\ — n 

V If - K a K bc 



(10) 



Here the semiclassical steady-state solution for the field 
amplitudes shows an abrupt transition from zero to 



T)J K c (Kb + Kbc) i 7 

a = - 2 Vf - 1 • 
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Po = \j 2 KaKc - Ve - 1 , 

V K a K bc 
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FIG. 3: (a) Mean occupation number in the photon laser (+) 
and the atom laser (o) above threshold. The lines illustrate 
the semiclassical solution while the inset shows the behavior 
of the atom laser around threshold, (b) The corresponding 
Mandel-Q parameter indicates sub-Poissonian statistics. 
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FIG. 4: Second-order correlation function g^ij) of the pho- 
ton laser (left) and the atom laser (right). 



e = p/pth, while 70 = p/k c below the threshold and 
remains constant at pth/K c above. Let n a and n b denote 
the mean occupation numbers of the two laser modes. 
When the atomic heating rate Kb c vanishes, we have 



K a KbK c 



(e — 1) , k = a,b, 



(12) 



such that both output laser beams have the same mean 
intensity I out = 2K a n a = 2K(,n&. As there is no other 
loss channel, this is even true for unequal decay rates K a 
and Kb- Only the heating rate Kbc determines the ra- 
tio of the output-field mean intensities and it holds that 
jout > jout Thig ig illustrate d m Fig. |3] (a) where n a and 

rib are plotted against the pumping strength /i. We chose 
the two output coupling rates to be equal (to k), so that 
the output intensities are related to the mode occupation 
numbers by the same scaling parameter. The other pa- 
rameters are (nbc, K c ,rj) = (0.1,10, 5)k and pth ~ 2.1k. 
The markers correspond to the solution of the master 
equation (J6j) obtained by a Monte Carlo wave function 
simulation. Each of the laser modes clearly exhibits an 
above-threshold population where the two solutions are 
in good agreement. 

In addition to average intensities we will now consider 
intensity fluctuations. For lasing the fluctuations of the 
intensity should be small compared to the mean. To show 
this we have plotted the Mandel-Q parameter, defined by 
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FIG. 5: (a) Noise spectrum V(fi) without heating (broken 
line) and with finite heating rate Kbc = k (full line), (b) 
Zero-frequency component of the noise spectrum. 

Q = (An 2 ) / (n) — 1, for the photon and atom statistics in 
Fig. 03(b). While it should be close to zero for a laser, the 
negative Q here even exhibits sub-Poissonian statistics. 
In addition the requirement of second-order coherence 
can be recast in terms of Glauber's second-order time 
correlation function as 

| ff ( 2 )(r)-l|«l. (13) 

Figure 02 depicts the coherence function g^ 2 \r) for each 
laser mode where both lasers operate above threshold 
and n a = Kb- 

As mentioned our system more closely related to the 
non-degenerate parametric oscillator than a laser. For 
photons the most striking consequence of parametric am- 
plification is noise reduction of quadrature fluctuations 
and intensity correlations below the level of vacuum fluc- 
tuations [20] . 

Since the "twin" signal and idler photons are created 
simultaneously by parametric down conversion, the fluc- 
tuations in the intensity difference are exactly zero in a 
system without dissipation. When cavity decay (which 
is a random process affecting both modes independently) 
is included, this zero-noise feature can still be recovered 
in the output photocurrent difference, provided that all 
photons are detected and the detection time is sufficiently 
large. Let us see what this means here. For ki, c = 
analogous arguments will hold for our system, while the 
additional dissipation channel due to heating of ground- 
state atoms will break the perfect symmetry in the out- 
put fields and suppress noise reduction. 

We now analyze the fluctuations in the output intensity 
difference AI = I° ut — I£ nt . The noise spectrum, which 
is defined as 




- [(AJ(r)AI(O)) - (AI) 2 ] , (14) 



is plotted in Fig. 03 (a). As expected it shows perfect 
noise reduction at zero frequency for Kb c — (broken 
line) while there is still squeezing of about 48% below the 
shot noise level for K{, c — k (full line). The zero- frequency 
contribution V(0) can be found in Fig. 0(b) as a function 



of Kbc- Hence the systems realizes a quantum correlated 
atom-photon pair source. 

In summary we showed that merging a photon laser 
and an atom laser in a cavity-confined optical lattice 
bears good prospects for a quantum source with gen- 
uine properties. The stimulated emission driven Raman 
side-band cooling leads to lower temperatures than cav- 
ity cooling and suppresses the severe problems of sponta- 
neous emission and reabsorption. For sufficient pumping 
both quantum modes exhibit coherent oscillations with 
sub-Poissonian output statistics. While the laser out- 
puts are second-order coherent, the output intensity dif- 
ference shows noise compression very similar to paramet- 
ric amplification and could be used to construct heralded 
coherent atom beams. We acknowledge funding by the 
Austrian Science Fund projects P17790 and S1512. 
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